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ABSTRACT 

The young star AB Aurigae is surrounded by a complex combination of gas-rich and dust dominated 
structures. The inner disk which has not been studied previously at sufficient resolution and imaging 
dynamic range seems to contain very little gas inside a radius of least 130 astronomical units (AU) 
from the star. Using adaptive-optics coronagraphy and polarimetry we have imaged the dust in an 
annulus between 43 and 302 AU from the star, a region never seen before. An azimuthal gap in an 
annulus of dust at a radius of 102 AU, along with a clearing at closer radii inside this annulus, suggests 
the formation of at least one small body at an orbital distance of about 100 AU. This structure seems 
consistent with crude models of mean motion resonances, or accumulation of material at two of the 
Lagrange points relative to the putative object and the star. We also report a low significance detection 
of a point source in this outer annulus of dust. This source may be an overdensity in the disk due 
to dust accreting onto an unseen companion. An alternate interpretation suggests that the object's 
mass is between 5 and 37 times the mass of Jupiter. The results have implications for circumstellar 
disk dynamics and planet formation. 

Subject headings: instrumentation: adaptive optics — methods: data analysis — stars: individual 
(HD31293) techniques: image processing — stars: planetary systems 



1. INTRODUCTION 

The star AB Aurigae (hereafter AB Aur; distance 
d = 144 parsecs; visual magnitude, V = 7.04 m ; spec- 
tral type AOVe; mass M = 2.4 ± 0.2 M ; age 1 - 3 
Myr; also known as HP 312 9 3, HIP 22910, BD+30 741; 
Ivan den Ancker et all (|1997t ): iDeWarf et all (|2003h ) has 
been intensively studied and is widely perceived as pro- 
totypical of young, intermediate-mass stars, the so-called 
Herbig Ae stars. The star belongs to the Taurus- Auriga 
star-forming region, which is thought to be between 1 and 
3 Myr old, where the dispersion in age s is most likely due 
to observational or calibration errors ([White et al.lf2007L 
e.g.). Because of its brightness, age and proximity to 
Earth, this star presents an important opportunity for 
investigating star and perhaps planet formation. 

Indeed, the star exhibits a fascinating array of phenom- 
ena, including a reflection nebula extending thousands of 
astronomical units (AU; the Ear th-Sun distance) away 
from the star (jGradv et al.lll999h : a more compact, an- 
nular (100-750 AU) region mainly composed of gas as 
detected in 13 CO, 12 CO, C 18 0, HC O+, H 2 , and millime- 
ter wavelength continuum imaging ((Mannings fc S argent! 
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1997 : iThi et all l200l iHenning et aD 119981 : iPietu et all 
20051 e.g.); and an apparent inner region within about 



100 AU of the star which exhibits gas depletion, where 
the ratio of the mass in dust to that in gas is large 
(base d on the surface densities derived by (|Pietu et all 
l2005h ). It seems that the dusty inner disk has an outer 
radius at about 130 AU (jFukagawa et all [200l l. out- 
side of which the gas column density rises. Also, near 
infrared (NIR) imaging reveals a spiral structure out- 
side of the inne r disk (from 130 to a few hundred AU; 
Fukag awa et al.1 (|2004l )). Ancillary to this study, but of 
importance in understanding this system, the material 
also exhi bits a central hole of rad ius between 0.3 and 
0.6 AU (jMillan-Gabet et al.l f200lh . A bright spot has 
also been detected about 1-4 AU away from the star 
(|Millan-Gabet et al.1 [20061 ) with NIR interferometry. 

Interpretations of these observations are inconsistent 
with each other in the literature, with some authors de- 
scribing disk- like structures and others not. The system 
is clearly complicated, and it may have properties of both 
gas-rich and debris disks, based on qualitative consid- 
eration of the observations mentioned above. As such, 
it may be in a stage of evolution that is intermediate 
between these two states, and it presents an ideal ob- 
servational opportunity to study disk evolution and the 
formation of small bodies in such disks. This means that 
high resolution observations of the region between 1 and 
130 AU are particularly important. 

However, this region has been difficult to image in light 
scattered from the dust because of the extreme bright- 
ness difference between the star and the dust. We present 
the first images of polarized light from this region in the 
NIR with sufficient contrast (residual starlight after sup- 
pression of < 10 -5 at radii of 0.3 to 2.2 arcseconds) and 
resolution (14 AU) to resolve structure in the dust disk. 
This level of starlight suppression was achieved using po- 
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Fig. 1. — Images of the stars HD 107146 (a and b) and AB Aur (c and d) showing the 4.7 arcsecond field of view with N up and E 
left, (a) and (c) are the summed intensity I images shown with logarithmic scale spanning the minimum and maximum counts. Significant 
starlight exists throughout the field of view for both stars, (b) and (d) are the double-differenced P images (see text) also shown with 
logarithmic color scale. The coronagraphic occulting spot is slightly smaller (0.45 arcsecond diameter) than the rejected region of each 
image, indicated by the shaded circle 0.6 arcsecond in diameter in both images. Cardinal direction arrows are 1 arcsecond long and drawn 
to scale. Note that the maximum count level in (d) is about ten times that in (b). 

ity to stars, because starlight can be efficiently removed 
by forming a polarized light imag e (iKuhn et al l l200lt 
iPerrin et al.ll2004j : lApai et al.ll2004l ; iTamura et alj|2006j) . 
Indeed, direct light from the star is only minimally po- 
larized by small-angle (< 5 x 10~ 4 degrees) forward scat- 
tering due to the Earth's atmosphere or the interstellar 
medium. In contrast, light from the star that has scat- 
tered off of material in the near-star environment is im- 
parted with far more polarization due to the large angle 
of scattering (near 90° for a face-on disk). Thus by differ- 
encing images obtained in orthogonal polarization states 
(i.e. I + Q and I - Q, I + U and I - U, or I + V and 
I—V, using the standard Stokes vector notation) the star 
light is removed, and an image of polarized light, a "P 



larimetry in conjunction with high-order adaptive optics 
(AO; for spatial resolution) and advanced coronagraphy 
(for star light suppression) u nder t he aegis of the Lyot 
Project (iQppenheimer et all |200l 1200a Iffinklev et all 
I2007t ISivaramakrishnan et al.ll200lh . 



2. POLARIMETRY 

Scattered light from the gas and dust in the vicinity of 
a star should be linearly polarized with the electric field 
perpendicular to the plane of incidence as seen from the 
telescope (assuming a single scattering process). Thus, 
dust in close proximity to a star will, in general, induce 
a spatially variable polarization angle and fractional po- 
larization as seen in a polarized light image. This makes 
imaging polarimetry of stars a valuable tool for study- 
ing star formation. In fact, numerous imaging polarime- 
try observations of disks have been made, permitting the 
observation of extremely faint sources in close proxim- 



image," can be formed where P = ^Q 2 + U 2 + V 2 . The 
spatial variation of the P image is a sensitive measure of 
variations in the scattering properties of the near-stellar 
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Fig. 2. — P image of the AB Aur system as in Fig. Id, but with a median filter applied with a five pixel bounding box, matched to half 
the diffraction-limited resolution of the instrument. The image at upper right is a magnified portion centered on the location of a point 
source (see text). In this box, the color scale is adjusted to aid interpretation. The color scale is logarithmic. 



dust and gas and is insensitive to atmospheric speckle 
noise. A further discussion can be found in the appendix 
of this work. The series of images acquired (six sets of 
two simultaneous images, each with different polarimeter 
settings; see Appendix) are also added together to form 
a separate total intensity or "J" image. 

3. OBSERVATIONS AND DATA REDUCTION 

Over a period of one hour and thirty-eight minutes, 
commencing at 09:55 on 2006 December 12 UTC, we ac- 
quired ten polarimetric measurements in H-band (1.45- 
1.83/im) with the Lyot Project coronograph installed 
on the 3.63m AEOS telescope atop Haleakala, Maui. 
Each measurement is a series of six dual-beam polar- 
ized images (60-s exposures) with different liquid-crystal 
retarder settings (see appendix for details), to permit 
derivation of the P image. The weather at the observa- 
tory was photometric, with winds of less than 8 km/h 
and relative humidity below 12%. The local Fried pa- 
rameter, tq, a measure of the strength of turbulence 



in the atmosphere above the observatory, spanned the 
range of 15 to 32 cm, indicating ideal conditions for 
AO observations at AEOS. AEOS has a 941-actuator 
AO system whic h suppresses imaging disto rtions due to 
the atmo sphere dRoberts &: Nevmanl I2002D. The coro- 
nagra ph (|Qppenheimer et al.ll2004l T2003; Hi nklev et al.l 
120071 ) is coupled with a NIR camera called "The Ker- 
mit" (Pe rrin et al.ll200"3l ), with which we image the 4.7 
arcsecond field of view centered on the star at a pixel 
scale of 13.7 mas pixel -1 . The diffraction-limited imag- 
ing resolution of the system is approximately 94 mas in 
.H-band. 

Used here as an experimental "control" observation, 
data on HD 107146 (a V = 7.07 m G2V star) were taken 
in a manner identical to that of AB Aur commencing at 
15:14 on 2006 December 12 UTC, with a similar range of 
weather conditions and telescope pointing relative to the 
horizon. We observed this star because it has a known 
debris disk, although the likelihood of our detecting it 
was extremely low, because the disk is optically very thin 
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and has only been imaged out side of the Lyot Project 
field of view (|Ardila et al.ll2004l ). Furthermore, a lack of 
infrared excess at 18 /im indicates a cent ral hole at least 
as w ide as 31 AU (1 arcsec) in radius (jWilliams et al.l 
12004( 1 . As expected, HD 107146's disk was not detected 
within our field of view (within 67 AU). Further discus- 
sion of these points can be found in the appendix. 

Calibration of the system includes standard dark 
current and flat-fielding gain images (acquired us- 
ing incandescent lamps each night), ba d pixel maps 
([Soummer et al.ll2006l; iHinklev et al.ll2007l ). and observa- 
tions of binary stars with well-known orbits to ascertain 
pixel scale and image rotation fiducials. The appendix 
to this work describes the steps in more detail. 

Fig. 1 shows the I and P images of both HD 107146 
and AB Aur. The structure we see in these images is 
better analyzed after filtering the image with a median 
pixel selection over a box matched to half the instru- 
ment's diffraction limit (3.5 pixels, rounded up to 5). 
Oversampling in the image plane aids speckle removal in 
postprocessing, but it can dilute real signals in final anal- 
yses. The median filtered P image of AB Aur is shown 
in Fig. 2. Before analyzing the structure in this image, 
we must establish quantitative sensitivity limits. 

4. DETECTION LIMITS 

According to a t echnique we developed previously 
(jHinklev et alJl2007l> . we have derived the difference in 
magnitude between the occulted star and a 3-cr point 
source as a function of position in the / images for HD 
107146, to evaluate the sensitivity of the AB Aur data 
without the presence of its disk, or any other source in 
the field of view. 

Fig. 3, a so - called "dynamic range plot" 
([Oppenheimer et al.l [2000), shows the difference in 
magnitude (AM) between the occulted star and a point 
source that would be detected at the 3-cr level in our I 
image, as a function of radial distance from the star. 
The photometric values were calibrated to unocculted 
images of AB Aur taken immediately prior to the 
occulted images. An extension of this technique was 
applied to the polarimetric data to derive the sensitivity 
limits for various fractional polarizations. For a real 
source in our data, only a lower limit to the fractional 
polarization is measurable, because of the contamination 
from the starlight throughout the image (see appendix). 

5. AB AUR IMAGES 

The images of the AB Aur system (Figs, lc, Id and 
2) exhibit several important features. In the I image, 
starlight dominates the entire field of view (Fig. lc) and 
shows some elongation along a position angle (PA) of 321 
±1°. This elongation, which is purely instrumental, and 
apparent in both the AB Aur data and the HD 107146 I - 
image, is due to a co mbination of the win d's effect on AO 
image improvement (|Sheehv et al.l [20061 . e.g.), and some 
residual image motion compensation errors. (The wind 
had an average prevailing direction of 320° during these 
observations.) In the P image (Fig. Id), light is detected 
throughout the field of view but is mostly constrained 
to a region 2.1 arcseconds (302 AU) in diameter centered 
near the star. Outside of the central region, the detection 
is of low significance, but does exhibit some amplification 
to the NE, which is where the brightest "spiral" feature 
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Fig. 3. — Dynamic Range Plot. This plot shows the difference in 
magnitude (AM) between a star occulted by the coronagraph and 
a 3-cr point source or surface brightness as a function of angular 
separation. Corresponding surface brightness is given on the right. 
The size of the occulting mask is indicated by the dark shaded 
region and the lighter shaded region corresponds to the rejected 
part of the data. The solid line indicates the dynamic range for 
the /-image while the dashed lines are dynamic range functions 
indicating the total intensity for sources with various fractional 
(%) polarizations. 

has been seen in other studies ([Fukagawa et al.ll2004D . In 
the inner region significant structure is present. The pho- 
tocenter of the inner disk is shifted 88 mas (12 AU) along 
a PA of 336 ± 2°. The outer annulus of material in this 
inner disk shows a distinctly non-uniform distribution in 
the azimuthal direction, with a deep depletion of polar- 
ization centered 700 mas (102 AU) away from the star at 
a PA of 3 33 ±2°. In a study c onducted at the Subaru 
telescope (Fukaga wa et al.l 12004') . there appears to be a 
similar depleted region in this outer annulus, although it 
is right on the edge of their coronagraphic mask. 

It is important to note that none of the structure seen 
in the P image correlates with structure in the I im- 
age, reinforcing the fact that only lower limits to the 
fractional P can be determined. For example, the I- 
image elongation is at least 15 degrees away from the 
depleted region's azimuthal angle in the outer annulus. 
The residual P/I amplitude in HD 107146 was about 
1%. In contrast the P/I structure we observe in AB Aur 
was between 1 and 9% — significantly larger than the 
residual systematic polarization noise we measure in the 
non-detection around HD 107146. 

In the center of the depleted region, there appears to 
be a 2.8-er amplification of polarization fraction spanning 
a seven-pixel diameter circular aperture ( ~ A/Z?Lyot, or 
94 mas). This appears to be an unresolved point source 
with P/I > 5%. We investigated all other features in the 
image and found none to have point source shape or scale 
(based on unocculted images) with significance greater 
than 2.5-er. No optical ghosts have ever been observed in 
the NIR in our system since it was commissioned in early 
2004, and we have never observed spurious reflections 
from bright stars outside the field of view. This putative 
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point source may or may not be real, and despite its 
marginal detection, we include it in our analysis below, 
although the principal observations of structure in the 
disk stand regardless of the reality of the point source. 
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Fig. 4. — Derived radial surface brightness of the AB Aur disk 
in Stokes /, using the HD 107146 data, where no polarized diffuse 
light was detected, as a reference azimuthally averaged brightness 
of an occulted point source. Error bars indicate one standard de- 
viation for the range of values observed in each annulus of data. 
The dashed line ind icates the radial profi le derived from Subaru 
coronagraphic data (Fukagawa ct al. 2 0041) . 

Further structure in the P image seems to indicate a 
drop in polarized light as one proceeds radially inward 
from the outer region toward the star. This is partic- 
ularly pronounced in the NW at a radius of about 570 
mas (83 AU). Closer to the star, the polarized light in- 
creases again up to the region obscured by the corona- 
graphic mask. There may be an annular depletion zone 
in close proximity to the coronagraphic mask edge, but 
we state this with less certainty. This annular depiction 
zone (which, also, shows no correlation with speckles in 
the I image) is located at roughly 300 mas (43 AU) ra- 
dius from the star and a PA of about 142°. 

6. INTERPRETATION: CIRCUMSTELLAR MATERIAL 

The azimuthally averaged surface brightness of the cir- 
cumstellar material provides initial measurements of the 
disk material. We use an azimuthal average, because 
the only effective way to remove the stellar light from 
the I-image, without restricting study to polarized light, 
is to use a reference star with no detectable circum- 
stellar material. AO point spread functions are highly 
variable and complex, meaning that point-spread func- 
tion subtraction is not practical (at best). An azimuthal 
average permits a first order understanding of the dif- 
ference between AB Aur and a reference star, and the 
simplest measure of flux from the disk alone. To this 
end, we used HD 107146's azimuthally averaged /-image 
subtracted from the same average for AB Aur (Fig. 4). 
The disk shows surface brightness ranging from about 80 
to 3 mJy/arcsec 2 at radii of 50 to 150 AU. This agrees 
with the values previously derived fro m 120 AU outward 
in the H-band (Fukaga wa et al.l [2004) . These values are 
also useable as a conservative upper limit to HD 107146's 
disk flux in the 8 to 66 AU region around that star. 

Interpretation of the morphology present in the P im- 
age depends on the notion that that polarization frac- 
tion in a given pixel in the image is a function of the 
column density of dust in the disk. What exactly this 



means in terms of the disk's structure is not entirely 
clear. A number of other studies suggest that this disk 
is o ptically thick in the m id-IR out to radii of about 118 
AU iMarifias et all (|2006f) . roughly the region we have 
detected, although data on the spatial scales reported 
here are lacking. If the disk is optically thick at the 
wavelengths we are observing here, then the P image re- 
veals only a map of the inhomogeneities in the surface of 
the disk, not the mid-plane features. If the disk is opti- 
cally thin at these wavelengths and spatial scales, then 
the P image will map the full structure of the dust in 
the disk. This is crit i cal to interpretation, because, as 
Uang-Condell fc Boss! ((2007) note, the surface structure 
in a relatively thick disk will not necessarily reflect the 
interior density structure. The full structure of the den- 
sity is necessary to utilize detailed modeling of disks to 
provide rigorous theoretical interpretation, therefore, we 
need to know whether the disk is optically thick or thin. 
However, the estimates of optical depth derived thus far 
are based on extremely simplistic models that assume az- 
imuthal and vertical homogeneity in the disk structure, 
with a radial dependence. What is beyond doubt is that 
these models are inaccurate based on the highly asym- 
metric morphology of the disk as revealed here. The 
conclusion must be that some parts of the disk may be 
optically thick while others may not be, including some 
of the regions imaged here. Numerical models, not sim- 
ple analytical ones are needed. 

The primary feature of the observed region is in the an- 
nulus at an approximate radius of 102 AU with a width 
equivalent to the image resolution ± 10 AU, where there 
is an apparent hole at a PA of 333°. In addition, there are 
clumps approximately 60° from the gap in the azimuthal 
direction, and a general reduction of density on the op- 
posite side of the gap. We note that the very outer edge 
of this depletion zone is seen in data reporting a spiral 
structure in the near-IR, right at the edge of the cor - 
responding coronagraphic mask ( Fukaga wa et al.l [20041. 
This structure, as noted above, also seems offset from the 
stellar position by about 12 AU. Since we are imaging in 
polarized light, some of this geometry could be due to 
the efficiency of scattering due to the inclination angle 
of the disk, variously reported for this inner region to 
be b etween 12° and 30° , with a major axis near 60 to 
80° ((M arihas et alj [20061. Indeed, if these numbers are 
correct, the annular depletion zone is nearly aligned with 
the minor axis of the disk and the opposite side would be 
expected to have reduced polarization fraction because it 
would be subject to the more weakly polarizing forward 
scattering by the dust. Because of the great uncertainty 
in the inclination angle and the large range in values 
depending on the material imaged in a given study, we 
did not see the value in adopting a particular inclination 
angle for the region we have imaged. Thus we have no 
"deprojected" image of the disk, which would also re- 
quire a model for the scattering material and its Mueller 
matrix. 

A veritable zoo of models for structure in circumstel- 
lar disks exists, with a dominant competition currently 
between so-called "core accretion" and "gravitational in- 
stability" theories to explain the formation of low mass 
objects in disks. These competing models predict differ- 
ent types of structure in disks. However, as far as could 
be determined from published material, the structure we 
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see in this disk is most similar to models where a small 
object is present tha t dynamically ind uces mean motion 
resonance structure (jWolf et all 120071 e.g). That struc- 
ture creates amplifications of density in the azimuthal di- 
rection roughly ±60° away from the for ming object, with 
a weaker clump opposite the object (jWolf et al.l I2007L 
esp. Fig. 10). This structure also has clearing just in- 
side the annulus as well as amplified density further in, 
another set of features that seem to exist in our P im- 
ages. There are issues with interpreting AB Aur's disk 
in this manner. These models exhibiting mean motion 
resonance generally require optically thin disks with no 
gas present. We made the point above that the optical 
depth is uncertain, and other observat ions show that thi s 
inner region is highly depleted of gas (jPietu et all feOOS). 
If this interpretation is correct, a small body would be 
forming in the gap and should be centered therein. How- 
ever, the theoretical modeling effort for disks of this type 
is not sufficiently mature to provide a complete interpre- 
tation. The observations can most simply be interpreted 
as showing an amplification of material around the L4 
and L5 Trojan or Lagrange points relative to an object 
that is situated in the gap. We discuss the possible de- 
tection of that object below. 

Comparing the AB Aur structure with other disks is 
important. First we note the large size of this disk and 
associated other structures (> 120 AU). A few low-mass 
companion formation theories exclude the formation of 
objects at such large distances as 102 AU, due to a pre- 
sump tion of very low dust density (jJang-Condell fc Boss! 
l2007l e.g. and references therein). These models claim 
that very low mass companions at such distances could 
only be there through dynamical interactions after form- 
ing near the star, and therefore should have highly eccen- 
tric orbits. The core accretion models can form objects 
at such distan ces, but generally of extremely low mass 
(Baraffc 2005), and interestingly, the gravitational in- 
stability models can only form objects in the 50 to 100 
AU range because the disk is too hot closer to the star 
(jJang-Condell fc Boss! 12007). However, observations of 
other disks reveal significant structure at scales certainly 
as large as this, albeit for older systems. For example, the 
star Fomalhaut has an offset ring of material, with some 
models claiming that this ring is due to a Saturn mass 
object at an orbital distance of ~ 120 AU (|Quilledl2006l) 
with low eccentricity. In some respects one could specu- 
late that the AB Aur structure is an evolutionary precur- 
sor to that found around the older Fomalhaut, which is 
very similar in physical properties to AB Aur except for 
its age. Furthermore, the star HD 142527 shows evidence 
for another offset (~ 20 AU) dis k of radius 500 AU an d 
a hole of about 125 AU in radius (Fukaga wa et al.| [2007). 
Such a structure is much larger than what we report here, 
and not too dissi milar to a smaller sc ale structure around 
epsilon Eridani (|Greaves et al.lll998h . 

Unfortunately few other Herbig Ae/Be stars have been 
observed at the spatial scales (< 100 AU) presented here 
or with polarimetry, so direct comparison with similarly 
aged objects is difficult. There are two possible excep- 
tions: observations of HD 141569 and HD 100546, which 
have been observed with, for example Hubble STIS on 
larger scales and in optic al scattered l i ght (which traces 
a variety of phenomena; iGradv et al.1 (|2005l 1200 ll )). In 
the case of HD 100546, complex structure is revealed at 



essentially all probed radii from the star and because it 
is significantly fainter at larger radii (> 500 AU) from 
the star than AB Aur, it may be more evolved and opti- 
cally thin. The observations of this star within 100 AU 
are, unfortunately lacking and so direct comparison with 
these observations provides little insight. 

For HD 1415 69, which is widely described as a "tran- 
sition al" disk (|Quillen et alj l2005t I Herrmann fc Krivovi 
I2007L e.g.), shows a spiral structure, also at large radii (> 
300 AU) and a continuous distribution o f dust at closer 
radii with multiple ring-like structures ( Clampin et al. 
2003). This has direct similarity with the structures we 
have observed in AB Aur, and the star is of similar spec- 
tra type and age. The azimuthal gap structure that we 
observe, though, does not seem to be observed in any 
other circumstellar disk material. 

On e other comparison can be drawn with the star 49 
Ceti (jWahhai et al.ll2007h . which is an A1V star about 
10-20 Myr old. Some authors claim that this object is 
evolving from the Herbig Ae stage to the more mature 
"Vega" state, and it is clear that there is a transition 
in dust properties with smaller grains dominating in the 
30-60 AU region dWahhai et al.l[200l . This may be a 
slightly more mature version of the disk we observed here 
and shows a similar radial hierarchy of structures and 
dust properties, even though in AB Aur the transition in 
dust properties seems to be at 130 AU, rather than at 60 
AU. The difference in age suggests that these very similar 
stars, when compared, provide clues to the evolution of 
disks around intermediate mass stars. 

7. INTERPRETATION: A PROBABLE COMPANION 

It may be possible to derive a mass for whatever object 
is causing the perturbations in the dust density in the 
disk through detailed modeling of the resonance struc- 
ture and an estimate of the mass of material in this struc- 
ture. This is beyond the scope of this paper, and requires 
significant theoretical work. Here, we explore the possi- 
bility that the point source reported above in the annular 
depletion zone is a bona fide companion of AB Aur, and 
provide two possible interpretations of what it is. This 
discussion assumes that the detection is real, even though 
it has low signal-to-noise ratio. 

Although we have not confirmed that the object de- 
tected in the gap is physically associated with AB Aur, 
through the detection of orbital motion or common 
proper motion, we have determined the probability of for- 
tuitous alignment with unassociated background sources. 
More details on this calculation are provided in the ap- 
pendix. The object's apparent H magnitude, m# , ranges 
from 12.76 in the brightest possibility where it is just be- 
low the detection limit in the I-image (at a radius of 
0.7 arcsec), and 17.26 in the faintest case where it is 
100% polarized. These are the sum of the values from 
the dynamic range cu rve and ran = 5. 06™ for AB Aur 
(|Skrutskie et al.l |2006[ ) . Our model indicates that the 
probability that a background object happens to lie in 
this location is less than 3.4 x 10~ 7 for the brighter limit 
and 9.9 x 10 -5 for the fainter one. The probability is so 
minute that we assume physical companionship. 

There are two ways to interpret this companion. First, 
we can derive limits to its mass by assuming that we 
have detected thermal emission from a cooling body 
coeval with the star. The i?-band (A = 1.6/xm) ab- 
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solute magnitude based on the distance to AB Aur is 
7.0 m < M H < 11.5™. Using a solar luminosity of 
Mb = 3.35 m , the H band luminosity of the object is 
between 3.5 xl0~ 2 and 5.6 xlO~ 4 L©, the luminosity of 
the Sun. Given the 1 to 3 Myr age of the system, this 
companion has a mass of be tween 5 and 3 7 Mj based 
on theoretical cooling models [Baraffe (20051). Although 
these models are not ideal for such young ages and the 
results are dependent on the initial conditions for the 
models, there is a lack of alternative theory to constrain 
these physical parameters. 

A second interpretation may be more appropriate in 
this case, because the object exhibits polarization of 
greater than 5%. Although planets orbiting stars should 
exhibit polar ization, when viewed in reflected light, as 
high as 50% (jTomasko et al.lll978l e.g.) there is no pos- 
sibility, due to the brightness, that this source represents 
light from the star reflected off of a tiny planet at 100 
AU. In fact, if it were reflected light, the simple geome- 
try of the reflection would demand that the object have a 
radius of 0.1 AU, about 10 times larger than the primary 
star, if the albedo were 1. The object is in the midst of a 
rather dusty region, albeit depleted from neighboring re- 
gions. Perhaps our images merely detect the infall of dust 
onto an obscured, accreting companion, rather than the 
companion itself. Models of planet formation suggest, in 
some cases, vortices of material surrounding the sites of 
cores of new planets. Such regions would be unresolved 
at the imaging resolution available to the Lyot Project, 
for example. As such, they would appear to be unre- 
solved clumps of dust that could scatter starlight just as 
the main part of the disk would. Other models, where 
the disk is denser (optically thick even), predict bright 
spots and shadows at or slightly offset from the location 
of planet accretion due to a hole formed by the influx 
of material onto the forming body (| Jang-Condelll [2007L 
e.g.). No models exist that permit us to estimate a mass 
for this object based on the polarimetry presented here. 
However, if this is an accretion process, one would cer- 
tainly expect a high degree of temporal variability as the 
companion is bombarded at random moments by larger 
particles or planetessimals as the disk continues to thin 
out. 

8. CONCLUSION 

AB Aur and its environment are indisputably a com- 
plex system that provide clues to the physical processes 
of star and planet formation. Our study provides some 
insight into a new scale on par with our own solar system. 
The combination of observations of this system suggest 
that (a) , because of the 1- 4 AU gap b etween the star and 
disk jMillan-Gabet et al.ll200ll . 120061 ). the star is mostly 
formed and has completed its main accretion phase, al- 
though there is significant outflow still; (b) planet for- 
mation, if it will happen, has already begun, due to the 
pres ence of an appare ntly organized disk with depleted 
gas ([Pietu et al.ll2005f ). some of which may be adsorbed 
onto the dust itself; (c) the region imaged here should 
have a similar temperature to the giant planet region 
of our solar system, and as such these data could pro- 
vide constraints on giant planet formation mechanisms; 



(d) the remaining thin dust and gas at larger radii are 
either flowing into or out of the disk, as indicated by 
strong spiral features at larger distances from the star 
tiFukagawa et al.ll2004t iGradv et al.lll999tlHenning et all 
119981 e.g.). In addition, projection effects complicate this 
interpretation particularly if there really is an organized 
inner disk (such as we have imaged), with external ma- 
terial distributed more isotropically, or at least not in a 
series of coplanar disks. 

Regarding the structure we have resolved, models 
of disks showing mean motion resonances also suggest 
strong outflows or infalls of material at coeval epochs 
with signifi cant material tra pped in mean motion reso- 
nant orbits (|Wolf et al.ll2007t ) . These are the only models 
that seem to come remotely close to mimicking the struc- 
ture we observe. We may be witnessing such a process in 
action for the first time, or some combination of processes 
not yet modeled, and the poorly detected point source 
in our data may be related to the structure in the disk. 
Higher resolution and more sensitive observations of this 
fascinating system may be able to constrain planet and 
star formation models and permit distinguishing between 
competing theory. Unfortunately such observations can 
only be expected once new, large sub-mm observatories, 
such as ALMA, are in operation. The critical need here 
involves spatial resolution superior to 100 mas and dy- 
namic range greater than 12 magnitudes. A major space 
mission such as the Terrestrial Planet Finder, would pro- 
vide a crystal clear view of what is happening in this fas- 
cinating system. In the near term, HST's NICMOS in- 
strument can conduct similar polarimetric observations 
at lower resolution, and possibly worse dynamic range. 
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APPENDIX 

ASTROMETRIC CALIBRATION AND DATA QUALITY 

We used the star HD 76943 for the astrometric calibration, observed at 14:13 on 2006 December 16 UTC under 
similar observing conditions. On this date, the orbit determination sho ws that the secondary component is at 0.546 
arcseconds at a position angle of 333.54° (N-E) (|Hartkopf k. Ma son 2000). Images of that star are shown in Fig. Al to 

demonstrate the fidelity of the AO correction, roughly estimated via the Strehl ratio, S = exp _<T , where a is the root 
mean square of the residual wave front error. S exc eeded 85% in these images, as determined through comparison 
with laboratory-based images with the coronagraph (|Roberts et al J 12004( 1 . This performance is similar in acquisition 
images of AB Aur taken immediately before the primary data acquisition. 




Fig. 5. — Image of the binary star HD 76943 used to determine plate scale and rotation parameters. The first Airy ring of the PSF 
is broken into four spots due to the oversized spider vanes in the coronagraph's Lyot stop. The secondary star has a separation of 0.546 
arcsec at PA 333.54°. Pixel scale is 13.7 mas and N is up with E to the left. 

POLARIMETRY 

iKuhn et al.l (|2001l ) describe a double -difference method that minimizes instrumental effects in imaging polarimetry. 
To derive our P images, following the IKuhn et al.l (|2001( 1 technique, we take differences of simultaneous orthogonal 
polarization images (e.g. I + Q and / — Q). We then difference these images again after spatially interchanging the 
polarization states IKuhn et alj (|2001f ) and then sum them in quadrature to form a "P" image P = \J Q 2 + U 2 + V 2 . 
Instrumental and other systematic polarization effects are removed or tend to be spatially constant in this polarized 
image. To achieve this, a Wollaston prism located immediately after the coronagraph's Lyot stop splits the light into 
two images, viewed simultaneously side-by-side on the Kermit's detector, with the two images carrying perpendicular 
vectors of the polarization. Using liquid crystal variable retarders (calibrated with a Glann- Taylor prism), the instru- 
ment can swap the vectors for each field of view or rotate them. As a result, we obtain six sets of two images. The two 
images are recorded on the same detector simultaneously and with the same exposure time. By swapping the vectors 
on the two images in a subsequent exposure — the two exposures comprising a "set" — we can remove instrumental 
effects induced by the different optical paths of the two simultaneous images. 

Taking the first two exposures, for example, we subtract the left and right images in both and then subtract the 
resulting final image. This provides an image in polarized light that we call "Q," with "U" and "V" images derived with 
the subsequent exposures in the same way. These Q, U and V images are from the reference frame of the instrument. 
In order to derive the Q, U, V images in the reference frame of the sky, they have to be reprojected taking into account 
the Mueller matrix of the telescope and instrument as a function of pointing in the sky. Howeve r, this Mueller matrix 
is not only extremely complex, but is also unknown (Ha rrington et all 120061 : [Harrin gton 2008). Fortunately, the P 
image is immune to these effects, because it is the quadrature sum of the Q, U and V images, but this "cross-talk" 
induced by the telescope requires that we also measure V, which is often excluded in observations similar to these. 

Our P image does not directly reveal the absolute polarization amplitude or direction of the near-star optical 
scattering, because the intensity image has significant spurious scattered starlight and atmospheric speckle noise 
(jHinklev et al.l [20071 . e.g.). In addition, as mentione d above, the AEOS telescope changes b oth the magnitude and 
polarization angle of incoming light to the telescope (|Harrington et aLl l2006: Ha rringtonl [20081 . Fortunately this spu- 
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rious polarization cross-talk is essentially constant over the telescope field-of-view, meaning that the relative structure 
in the P image is not instrumental in nature. The observed AB Aurigae Q, U and V images were all very similar 
in nature with small relative rotations in the polarization angle and degree of polarization as the telescope changes 
pointing. However, the quantity \/ Q 2 + U 2 + V 2 was constant during all ten polarization sequences (6 dual images 
for each). Despite these limitations, the spatial variation of the P image is a sensitive measure of variations in the 
scattering properties of the near-stellar dust and gas and is insensitive to atmospheric speckle noise. 

EFFICACY OF SPECKLE SUPPRESSION 

To demonstrate how completely the starlight (and speckle noise) is removed with our imaging polarimetry process, 
we examine, in more detail, the observations of the star HD 107146 which is used in the main article as a "control" 
non-detection. It is bright in the coronagraphic I image (Fig. la), with a total of 1.33 xlO 8 photons detected in a 360 
second sequence of exposures. However, the star is almost entirely absent in the P image (Fig. lb) with 1.60xl0 5 
photons, most of which, due to diffraction, are behind the coronagraphic occulting mask, a region of no consequence 
to scientific results. The polarimetry method reduces the residual starlight by a factor of 843 in this case and brings 
the vast majority of the field of view to the instrumental throughput limits — the P image shows a typical readnoise- 
limited background. No speckle structure remains in the region usable for scientific investigations. Combined with 
the coronagraph, the central star's light has been suppressed by a factor of about 10 -5 over the entire field of view, 
and possibly by much more since our images were not limited by the starlight but rather instrument sensitivity. It is 
important to note that none of the speckles in the / image are present in the P image. This is true for the AB Aur 
data as well. 

POLARIMETRY SENSITIVITY 

To determine our sensitivity in the P image, artificial point sources with a given intensity and fractional polarization 
were placed radially into the constituent images prior to composing the P image, with appropriate adjustments in the 
pairs of simultaneous images so that the relevant fractional polarization was correctly represented by the fake sources. 
These images were double differenced as described above (Polarimetry) . The total intensity of these point sources in 
the original images was adjusted until a 3-cr detection was obtained in the processed Q, U , and V data. (P images 
have non-Gaussian noise properties so noise must be derived from the constituent images.) The total intensity (not 
just the flux in net-polarized light) for a given fractional polarization is plotted in the dynamic range plot, although for 
any real source detected, only a lower limit to the fractional polarization is measurable, because of the contamination 
from the starlight throughout the image. 

PROBABILITY OF POINT SOURCE CONTAMINATION 

We constrained the probability that a point source would be co-aligned with features in AB Aur's disk using our 
detection limits and a model for the density of object s on the sky. Deriving a star count-magnitude relationship from 
the Two- Micron All-Sky Su rvey Point Source Catalog ()Skrutskie et al.ll2006f l and using a fit to galaxy count-magnitude 
relations (|Chen et al.l l2002). we find that sources in the region of sky around AB Aur have a sky surface density of 

8.73 + 0.30m H (El) 

4.69 + 0.45m H (E2) 

where N s and N g are the number of stars and galaxies per square degree, respectively, and mjj is the apparent H-band 
magnitude. 

Because we have detected this source in polarized light, we must modify these equations to reflect the f act that the 
majority of celestial sources have imme asurable (< 1%) po larization with our technique (|Whittetll2004h . Based on 
results of a survey of polarized starlight (jFosalba et al.ll2002h . we very conservatively assume that 50% of all stars and 
galaxies are polarized to the point where they would be detected in our P-image. This adds a constant of -0.30 to the 
equations above. 

To apply this, we make two estimates of the object's mjj, which are 12.76 in the brightest possibility where it is 
just below the detection limit in the /-image (at a radius of 0.7 arcsec), and 17.26 in the faintest case where it is 100% 
polar ized. These are the sum of the values from the dynamic range curve and m# = 5.06 m for AB Aur (|Skrutskie et al.l 
2006). From equations (1) and (2) and the area of the depleted region in the outer disk annulus, a circle of radius 0.2 
arcsec, the probability that a background object happens to lie in this location is less than 3.4 x 10~ 7 for the brighter 
limit and 9.9 x 10~ 5 for the fainter one. The probability is so minute that we assume physical companionship. 



log 
log 



dN a 
dm,H 



10 



Oppenheimer et al. 



REFERENCES 



Apai, D., Pascucci, I., Brandner, W., Henning, T., Lenzen, R., 
Potter, D. E., Lagrange, A.-M., & Rousset, G. 2004, Astron. 
Astrophys., 415, 671 

Ardila, D. R., Golimowski, D. A., Krist, J. E., Clampin, M., 
Williams, J. P., Blakeslee, J. P., Ford, H. C, Hartig, G. F., & 
Illingworth, G. D. 2004, Astrophys. J. Lett., 617, L147 

BarafTc, I. 2005, Space Science Reviews, 116, 67 

Chen, H.-W., McCarthy, P. J., Marzke, R. O., Wilson, J., Carlberg, 
R. G., Firth, A. E., Persson, S. E., Sabbey, C. N., Lewis, J. R., 
McMahon, R. G., Lahav, O., Ellis, R. S., Martini, P., Abraham, 
R. G., Ocmlcr, A., Murphy, D. C, Somerville, R. S., Beckett, 
M. G., & Mackay, C. D. 2002, Astrophys. J., 570, 54 

Clampin, M., Krist, J. E., Ardila, D. R., Golimowski, D. A., Hartig, 
G. F., Ford, H. C, Illingworth, G. D., Bartko, F., Bem'tez, N., 
Blakeslee, J. P., Bouwens, R. J., Broadhurst, T. J., Brown, R. A., 
Burrows, C. J., Cheng, E. S., Cross, N. J. G., Feldman, P. D., 
Franx, M., Gronwall, C, Infante, L., Kimble, R. A., Lesser, 
M. P., Martcl, A. R., Menanteau, F., Meurer, G. R., Milcy, 

G. K., Postman, M., Rosati, P., Sirianni, M., Sparks, W. B., 
Tran, H. D., Tsvetanov, Z. I., White, R. L., & Zheng, W. 2003, 
AJ, 126, 385 

DeWarf, L. E., Sepinsky, J. F., Guinan, E. F., Ribas, I., & Nadalin, 
I. 2003, Astrophys. J., 590, 357 

Fosalba, P., Lazarian, A., Prunet, S., & Tauber, J. A. 2002, 
Astrophys. J., 564, 762 

Fukagawa, M., Hayashi, M., Tamura, M., Itoh, Y., Hayashi, S. S., 
Oasa, Y., Takeuchi, T., Morino, J.-i., Murakawa, K., Oya, S., 
Yamashita, T., Suto, H., Mayama, S., Naoi, T., Ishii, M., Pyo, 
T.-S., Nishikawa, T., Takato, N., Usuda, T., Ando, H., Iye, M., 
Miyama, S. M., & Kaifu, N. 2004, Astrophys. J. Lett., 605, L53 

Fukagawa, M., Itoh, Y., Oasa, Y., Kudo, T., Fujiwara, H., Tamura, 
M., Hayashi, M., & Hayashi, S. 2007, in In the Spirit of Bernard 
Lyot: The Direct Detection of Planets and Circumstcllar Disks 
in the 21st Century, ed. P. Kalas 

Grady, C. A., Polomski, E. F., Henning, T., Stccklum, B., 
Woodgate, B. E., Telesco, C. M., Pina, R. K., Gull, T. R., 
Boggess, A., Bowers, C. W., Bruhweiler, F. C, Clampin, M., 
Danks, A. C, Green, R. F., Heap, S. R., Hutchings, J. B., 
Jenkins, E. B., Joseph, C, Kaiser, M. E., Kimble, R. A., 
Kraemer, S., Lindler, D., Linsky, J. L., Maran, S. P., Moos, 

H. W., Plait, P., Rocsler, F., Timothy, J. G., & Weistrop, D. 
2001, AJ, 122, 3396 

Grady, C. A., Woodgate, B., Bruhweiler, F. C, Boggess, A., Plait, 
P., Lindler, D. J., Clampin, M., & Kalas, P. 1999, Astrophys. J. 
Lett., 523, L151 

Grady, C. A., Woodgate, B. E., Bowers, C. W., Gull, T. R., Sitko, 
M. L., Carpenter, W. J., Lynch, D. K., Russell, R. W., Perry, 
R. B., Williger, G. M., Roberge, A., Bouret, J.-C, & Sahu, M. 
2005, ApJ, 630, 958 

Greaves, J. S., Holland, W. S., Moriarty-Schieven, G., Jenness, T., 
Dent, W. R. F., Zuckerman, B., McCarthy, C, Webb, R. A., 
Butner, H. M., Gear, W. K., & Walker, H. J. 1998, Astrophys. 
J., 506, L133 

Harrington, D. M., Kuhn, J. R., & Whitman, K. 2006, Publ. 

Astron. Soc. Pac, 118, 845 
Harrington, D. M. & Kuhn, J. R. 2008, PASP, 120, in press 
Hartkopf, W. I. & Mason, B. D. 2000, in Bulletin of the American 

Astronomical Society, 870 — h 
Henning, T., Burkcrt, A., Launhardt, R., Leinert, C, & Stccklum, 

B. 1998, Astron. Astrophys., 336, 565 
Herrmann, F. & Krivov, A. V. 2007, A&A, 476, 829 
Hinklcy, S., Oppenheimer, B. R., Soummer, R., Sivaramakrishnan, 

A., Roberts, Jr., L. C, Kuhn, J., Makidon, R. B., Perrin, M. D., 

Lloyd, J. P., Kratter, K., & Brenner, D. 2007, Astrophys. J., 654, 

633 

Jang-Condcll, H. 2007, in Proceedings of the conference In the 
Spirit of Bernard Lyot: The Direct Detection of Planets and 
Circumstellar Disks in the 21st Century. June 04 - 08, 2007. 
University of California, Berkeley, CA, USA. Edited by Paul 
Kalas., ed. P. Kalas 
Jang-Condell, H. & Boss, A. P. 2007, Astrophys. J. Lett., 659, L169 
Kuhn, J. R., Potter, D., & Parise, B. 2001, Astrophys. J. Lett., 
553, L189 

Mannings, V. & Sargent, A. I. 1997, Astrophys. J., 490, 792 
Marinas, N., Telesco, C. M., Fisher, R. S., Packham, C, & 
Radomski, J. T. 2006, Astrophys. J., 653, 1353 



Millan-Gabet, R., Monnier, J. D., Berger, J. -P., Traub, W. A., 
Schloerb, F. P., Pedretti, E., Benisty, M., Carleton, N. P., 
Hagucnaucr, P., Kern, P., Labcye, P., Lacasse, M. G., Malbet, 

F. , Perraut, K., Pearlman, M., & Thureau, N. 2006, Astrophys. 
J. Lett., 645, L77 

Millan-Gabet, R., Schloerb, F. P., & Traub, W. A. 2001, Astrophys. 
J., 546, 358 

Oppenheimer, B. R., Dekany, R. G., Hayward, T. L., Brandl, B., 
Troy, M., & Bloemhof, E. E. 2000, in Proc. SPIE, Vol. 4007, 
Adaptive Optical Systems Technology, ed. P. L. Wizinowich, 
899-905 

Oppenheimer, B. R., Digby, A. P., Newburgh, L., Brenner, 
D., Shara, M., Mey, J., Mandeville, C, Makidon, R. B., 
Sivaramakrishnan, A., Soummer, R., Graham, J. R., Kalas, P., 
Perrin, M. D., Roberts, Jr., L. C, Kuhn, J. R., Whitman, K., 
& Lloyd, J. P. 2004, in Proc. SPIE, Vol. 5490, Advancements 
in Adaptive Optics, ed. D. Bonaccini Calia, B. L. Ellerbroek, & 
R. Ragazzoni, 433-442 

Oppenheimer, B. R., Sivaramakrishnan, A., & Makidon, R. B. 
2003, Imaging Exoplanets: The Role of Small Telescopes, Vol. 
Ill (The Future of Small Telescopes In The New Millennium), 
155-181 

Perrin, M. D., Graham, J. R., Kalas, P., Lloyd, J. P., Max, C. E., 
Gavel, D. T., Pennington, D. M., & Gates, E. L. 2004, Science, 
303, 1345 

Perrin, M. D., Graham, J. R., Trumpis, M., , Kuhn, J. R., 
Whitman, K., Coulter, R., , Lloyd, J. P., & Roberts, L. C. 

2003, in 2002 AMOS Technical Conference, P. W. Kcrvin, J. 
L. Africano; cds. 

Pietu, V., Guillotcau, S., & Dutrey, A. 2005, Astron. Astrophys., 
443, 945 

Quillen, A. C. 2006, Mon. Not. Royal Astron. Soc., 372, L14 
Quillen, A. C, Varniere, P., Minchev, I., & Frank, A. 2005, AJ, 
129, 2481 

Roberts, Jr., L. C. & Neyman, C. R. 2002, Publ. Astron. Soc. Pac., 
114, 1260 

Roberts, Jr., L. C, Perrin, M. D., Marchis, F., Sivaramakrishnan, 

A. , Makidon, R. B., Christou, J. C, Macintosh, B. A., Poynecr, 
L. A., van Dam, M. A., & Troy, M. 2004, in Proc. SPIE, Vol. 
5490, Advancements in Adaptive Optics, ed. D. Bonaccini Calia, 

B. L. Ellerbroek, & R. Ragazzoni, 504-515 

Sheehy, C. D., McCrady, N., & Graham, J. R. 2006, Astrophys. J., 
647, 1517 

Sivaramakrishnan, A., Korcsko, C. D., Makidon, R. B., Bcrkcfcld, 
T., & Kuchner, M. J. 2001, ApJ, 552, 397 

Skrutskic, M. F., Cutri, R. M., Stiening, R., Weinberg, M. D., 
Schneider, S., Carpenter, J. M., Beichman, C, Capps, R., 
Chester, T., Elias, J., Huchra, J., Liebert, J., Lonsdale, C, 
Monet, D. G., Price, S., Seitzer, P., Jarrett, T., Kirkpatrick, 
J. D., Gizis, J. E., Howard, E., Evans, T., Fowler, J., Fullmer, 
L., Hurt, R., Light, R., Kopan, E. L., Marsh, K. A., McCallon, 
H. L., Tarn, R., Van Dyk, S., & Wheelock, S. 2006, Astron. J., 
131, 1163 

Soummer, R., Oppenheimer, B. R., Hinkley, S., Sivaramakrishnan, 

A. , Makidon, R. B., Digby, A., Brenner, D., Kuhn, J., Perrin, 
M. D., Roberts, Jr., L. C, & Kratter, K. 2006, in EAS 
Publications Series, 199-212 

Tamura, M., Fukagawa, M., Kimura, H., Yamamoto, T., Suto, H., 

k. Abe, L. 2006, Astrophys. J., 641, 1172 
Thi, W. F., Blake, G. A., van Dishocck, E. F., van Zadclhoff, G. J., 

Horn, J. M. M., Becklin, E. E., Mannings, V., Sargent, A. I., van 

den Ancker, M. E., & Natta, A. 2001, Nature, 409, 60 
Tomasko, M. G., West, R. A., & Castillo, N. D. 1978, Icarus, 33, 

558 

van den Ancker, M. E., The, P. S., Tjin A Djie, H. R. E., Catala, 

C. de Winter, D., Blondel, P. F. C, & Waters, L. B. F. M. 1997, 
Astron. Astrophys., 324, L33 

Wahhaj, Z., Koerner, D. W., & Sargent, A. I. 2007, ApJ, 661, 368 
White, R. J., Greene, T. P., Doppmann, G. W., Covey, K. R., 
& Hillenbrand, L. A. 2007, in Protostars and Planets V, ed. 

B. Rcipurth, D. Jewitt, & K. Kcil, 117-132 

Whittet, D. C. B. 2004, in Astronomical Society of the Pacific 
Conference Series, Vol. 309, Astrophysics of Dust, ed. A. N. Witt, 

G. C. Clayton, & B. T. Draine, 65-+ 

Williams, J. P., Najita, J., Liu, M. C, Bottinelli, S., Carpenter, 
J. M., Hillenbrand, L. A., Meyer, M. R., & Soderblom, D. R. 

2004, Astrophys. J., 604, 414 



The Solar-System-Scale Disk Around AB Aurigae 



Wolf, S., Moro-Martm, A., & D'Angelo, G. 2007, Planetary and 
Space Science, 55, 569 



